Raw data for "Linkage disequilibrium estimation of effective population size with immigrants from divergent populations: a case study on Spanish mackerel (Scomberomorus commerson)."

Abstract of Macbeth, G. M., Broderick, D., Buckworth, R. & Ovenden, J. R. (In press, Feb 2013). Linkage disequilibrium estimation of effective population size with immigrants from divergent populations: a case study on Spanish mackerel (Scomberomorus commerson). G3: Genes, Genomes and Genetics. Estimates of genetic effective population size (Ne) using molecular markers are a potentially useful tool for the management of endangered through to commercial species. But, pitfalls are predicted when the effective size is large, as estimates require large numbers of samples from wild populations for statistical validity. Our simulations showed that linkage disequilibrium estimates of Ne up to 10,000 with finite confidence limits can be achieved with sample sizes around 5000. This was deduced from empirical allele frequencies of seven polymorphic microsatellite loci in a commercially harvested fisheries species, the narrow barred Spanish mackerel (Scomberomorus commerson). As expected, the smallest standard deviation of Ne estimates occurred when low frequency alleles were excluded. Additional simulations indicated that the linkage disequilibrium method was sensitive to small numbers of genotypes from cryptic species or conspecific immigrants. A correspondence analysis algorithm was developed to detect and remove outlier genotypes that could possibly be inadvertently sampled from cryptic species or non-breeding immigrants from genetically separate populations. Simulations demonstrated the value of this approach in Spanish mackerel data. When putative immigrants were removed from the empirical data, 95% of the Ne estimates from jacknife resampling were above 24,000.

Genetic variation within two sympatric spotted gum eucalypts exceeds between taxa variation

Population substructure and hybridization, among other factors, have the potential to cause erroneous associations in linkage disequilibrium (LD) mapping. Two closely related spotted gum eucalypts, Corymbia variegata and C. henryi (Myrtaceae) occur in sympatry in the east coast of Australia and potentially interbreed. They are morphologically similar but are distinguished as separate species based on capsule and foliage size. To determine whether they hybridize in nature and its implications for LD mapping, we investigated the level of molecular divergence between the two species at two sympatric locations separated by 300 kilometres. Very few individuals of intermediate morphology were identified, despite the two species occurring only metres apart. Analysis of genetic structure using 12 microsatellite loci showed that genetic differentiation between populations of the same species at different locations (FST = 0.07 for both species; p = 0.0001) was significantly higher than that observed between species at each location (mean FST = 0.02 and 0.04 for Cherry tree and Bunyaville respectively; p = 0.0001; all Mann-Whitney U-test p ≤ 0.01). No species-specific alleles or significant allele frequency differences were detected within a site, suggesting recurrent local gene flow between the two species. The lack of significant allele frequency differences implies no population stratification along taxonomic lines. This suggested that there is little concern for cryptic hybridization when sampling from sites of sympatry for LD mapping.

Quantitative and molecular genetics of juvenile wood traits in radiata and slash/Caribbean pines.

The Juvenile Wood Initiative (JWI) project has been running successfully since July 2003 under a Research Agreement with FWPA and Letters of Association with the consortium partners STBA (Southern Tree Breeding Association), ArborGen and FPQ (Forestry Plantations Queensland). Over the last five and half years, JWI scientists in CSIRO, FPQ, and STBA have completed all 12 major milestones and 28 component milestones according to the project schedule. We have made benchmark progress in understanding the genetic control of wood formation and interrelationships among wood traits. The project has made 15 primary scientific findings and several results have been adopted by industry as summarized below. This progress was detailed in 10 technical reports to funding organizations and industry clients. Team scientists produced 16 scientific manuscripts (8 published, 1 in press, 2 submitted, and several others in the process of submission) and 15 conference papers or presentations. Primary Scientific Findings. The 15 major scientific findings related to wood science, inheritance and the genetic basis of juvenile wood traits are: 1. An optimal method to predict stiffness of standing trees in slash/Caribbean pine is to combine gravimetric basic density from 12 mm increment cores with a standing tree prediction of MoE using a time of flight acoustic tool. This was the most accurate and cheapest way to rank trees for breeding selection for slash/Caribbean hybrid pine. This method was also recommended for radiata pine. 2. Wood density breeding values were predicted for the first time in the STBA breeding population using a large sample of 7,078 trees (increment cores) and it was estimated that selection of the best 250 trees for deployment will produce wood density gains of 12.4%. 3. Large genetic variation for a suite of wood quality traits including density, MFA, spiral grain, shrinkage, acoustic and non-acoustic stiffness (MoE) for clear wood and standing trees were observed. Genetic gains of between 8 and 49% were predicted for these wood quality traits with selection intensity between 1 to 10% for radiata pine. 4. Site had a major effect on juvenile-mature wood transition age and the effect of selective breeding for a shorter juvenile wood formation phase was only moderate (about 10% genetic gain with 10% selection intensity, equivalent to about 2 years reduction of juvenile wood). 5. The study found no usable site by genotype interactions for the wood quality traits of density, MFA and MoE for both radiata and slash/Caribbean pines, suggesting that assessment of wood properties on one or two sites will provide reliable estimates of the genetic worth of individuals for use in future breeding. 6. There were significant and sizable genotype by environment interactions between the mainland and Tasmanian regions and within Tasmania for DBH and branch size. 7. Strong genetic correlations between rings for density, MFA and MoE for both radiata and slash/Caribbean pines were observed. This suggests that selection for improved wood properties in the innermost rings would also result in improvement of wood properties in the subsequent rings, as well as improved average performance of the entire core. 8. Strong genetic correlations between pure species and hybrid performance for each of the wood quality traits were observed in the hybrid pines. Parental performance can be used to identify the hybrid families which are most likely to have superior juvenile wood properties of the slash/Caribbean F1 hybrid in southeast Queensland. 9. Large unfavourable genetic correlations between growth and wood quality traits were a prominent feature in radiata pine, indicating that overcoming this unfavourable genetic correlation will be a major technical issue in progressing radiata pine breeding. 10. The project created the first radiata pine 18 k cDNA microarray and generated 5,952 radiata pine xylogenesis expressed sequence tags (ESTs) which assembled into 3,304 unigenes. 11. A total of 348 genes were identified as preferentially expressed genes in earlywood or latewood while a total of 168 genes were identified as preferentially expressed genes in either juvenile or mature wood. 12. Juvenile earlywood has a distinct transcriptome relative to other stages of wood development. 13. Discovered rapid decay of linkage disequilibrium (LD) in radiata pine with LD decaying to approximately 50% within 1,700 base pairs (within a typical gene). A total of 913 SNPS from sequencing 177,380 base pairs were identified for association genetic studies. 14. 149 SNPs from 44 genes and 255 SNPs from a further 51 genes (total 95 genes) were selected for association analysis with 62 wood traits, and 30 SNPs were shortlisted for their significant association with variation of wood quality traits (density, MFA and MoE) with individual significant SNPs accounting for between 1.9 and 9.7% of the total genetic variation in traits. 15. Index selection using breeding objectives was the most profitable selection method for radiata pine, but in the long term it may not be the most effective in dealing with negative genetic correlations between wood volume and quality traits. A combination of economic and biological approaches may be needed to deal with the strong adverse correlation.

Fifteen microsatellite loci for the jungle perch, Kuhlia rupestris

We developed and optimized 15 polymorphic microsatellite loci in the jungle perch, Kuhlia rupestris. Loci were screened in a single population (n = 24) from Fraser Island, Queensland, Australia. Number of alleles per locus ranged from 3 to 19 and observed heterozygosity from 0.25 to 1. No significant linkage disequilibrium was detected between any pair of loci. Genotype proportions for these loci in the population sampled were in Hardy–Weinberg equilibrium.

Genetic Structure of Cochliobolus sativus Populations Sampled from Roots and Leaves of Barley and Wheat in North Dakota

Common root rot (CRR) and spot blotch, caused by Cochliobolus sativus (Ito and Kurib.) Drechsl. ex Dast., are important diseases of barley (Hordeum vulgare L.) and wheat (Triticum aestivum L.) worldwide. However, the population biology of C. sativus is still poorly understood. In this study, the genetic structure of three C. sativus populations, consisting of isolates sampled respectively from barley leaves (BL), barley roots (BR) and wheat roots (WR) in North Dakota, was analysed with amplified fragment length polymorphism (AFLP) markers. A total of 127 AFLP loci were generated among 208 C. sativus isolates analysed with three primer combinations. Gene diversity (H = 0.277-0.335) were high in all three populations. Genetic variation among C. sativus individuals within population accounted for 74%, whereas 26% of the genetic variation was explained among populations. Genetic differentiation was high (empty set PT = 0.261, corrected G ''(st)= 0.39), whereas gene flow (Nm) ranged from 1.27 to 1.56 among the three populations analysed. The multilocus linkage disequilibrium (LD) ((r) over bard = 0.0760.117) was moderate in C. sativus populations. Cluster analyses indicate that C. sativus populations differentiated according to the hosts (barley and wheat) and tissues (root and leaf) although generalists also exist in North Dakota. Crop breeding may benefit from combining genes for resistance against both specialists and generalists of C. sativus.

Association genetics in Corymbia citriodora subsp variegata identifies single nucleotide polymorphisms affecting wood growth and cellulosic pulp yield

Wood is an important biological resource which contributes to nutrient and hydrology cycles through ecosystems, and provides structural support at the plant level. Thousands of genes are involved in wood development, yet their effects on phenotype are not well understood. We have exploited the low genomic linkage disequilibrium (LD) and abundant phenotypic variation of forest trees to explore allelic diversity underlying wood traits in an association study. Candidate gene allelic diversity was modelled against quantitative variation to identify SNPs influencing wood properties, growth and disease resistance across three populations of Corymbia citriodora subsp. variegata, a forest tree of eastern Australia. Nine single nucleotide polymorphism (SNP) associations from six genes were identified in a discovery population (833 individuals). Associations were subsequently tested in two smaller populations (130160 individuals), validating our findings in three cases for actin 7 (ACT7) and COP1 interacting protein 7 (CIP7). The results imply a functional role for these genes in mediating wood chemical composition and growth, respectively. A flip in the effect of ACT7 on pulp yield between populations suggests gene by environment interactions are at play. Existing evidence of gene function lends strength to the observed associations, and in the case of CIP7 supports a role in cortical photosynthesis.

Genetic Structure of Cochliobolus sativus Populations Sampled from Roots and Leaves of Barley and Wheat in North Dakota

Common root rot (CRR) and spot blotch, caused by Cochliobolus sativus (Ito and Kurib.) Drechsl. ex Dast., are important diseases of barley (Hordeum vulgare L.) and wheat (Triticum aestivum L.) worldwide. However, the population biology of C. sativus is still poorly understood. In this study, the genetic structure of three C. sativus populations, consisting of isolates sampled respectively from barley leaves (BL), barley roots (BR) and wheat roots (WR) in North Dakota, was analysed with amplified fragment length polymorphism (AFLP) markers. A total of 127 AFLP loci were generated among 208 C. sativus isolates analysed with three primer combinations. Gene diversity (H = 0.277-0.335) were high in all three populations. Genetic variation among C. sativus individuals within population accounted for 74%, whereas 26% of the genetic variation was explained among populations. Genetic differentiation was high (empty set PT = 0.261, corrected G ''(st)= 0.39), whereas gene flow (Nm) ranged from 1.27 to 1.56 among the three populations analysed. The multilocus linkage disequilibrium (LD) ((r) over bard = 0.0760.117) was moderate in C. sativus populations. Cluster analyses indicate that C. sativus populations differentiated according to the hosts (barley and wheat) and tissues (root and leaf) although generalists also exist in North Dakota. Crop breeding may benefit from combining genes for resistance against both specialists and generalists of C. sativus.

NEESTIMATOR v2: re-implementation of software for the estimation of contemporary effective population size (N-e) from genetic data

NeEstimator v2 is a completely revised and updated implementation of software that produces estimates of contemporary effective population size, using several different methods and a single input file. NeEstimator v2 includes three single-sample estimators (updated versions of the linkage disequilibrium and heterozygote-excess methods, and a new method based on molecular coancestry), as well as the two-sample (moment-based temporal) method. New features include the following: (i) an improved method for accounting for missing data; (ii) options for screening out rare alleles; (iii) confidence intervals for all methods; (iv) the ability to analyse data sets with large numbers of genetic markers (10000 or more); (v) options for batch processing large numbers of different data sets, which will facilitate cross-method comparisons using simulated data; and (vi) correction for temporal estimates when individuals sampled are not removed from the population (Plan I sampling). The user is given considerable control over input data and composition, and format of output files. The freely available software has a new JAVA interface and runs under MacOS, Linux and Windows.

Construction of microsatellite linkage maps for Corymbia

The genus Corymbia is closely related to the genus Eucalyptus, and like Eucalyptus contains tree species that are important for sub-tropical forestry. Corymbia's close relationship with Eucalyptus suggests genetic studies in Corymbia should benefit from transfer of genetic information from its more intensively studied relatives. Here we report a genetic map for Corymbia spp. based on microsatellite markers identified de novo in Corymbia sp or transferred from Eucalyptus. A framework consensus map was generated from an outbred F 2 population (n = 90) created by crossing two unrelated Corymbia torelliana x C. citriodora subsp. variegata F1 trees. The map had a total length of 367 cM (Kosambi) and was composed of 46 microsatellite markers distributed across 13 linkage groups (LOD 3). A high proportion of Eucalyptus microsatellites (90%) transferred to Corymbia. Comparative analysis between the Corymbia map and a published Eucalyptus map identified eight homeologous linkage groups in Corymbia with 13 markers mapping on one or both maps. Further comparative analysis was limited by low power to detect linkage due to low genome coverage in Corymbia, however, there was no convincing evidence for chromosomal structural differences because instances of non-synteny were associated with large distances on the Eucalyptus map. Segregation distortion was primarily restricted to a single linkage group and due to a deficit of hybrid genotypes, suggesting that hybrid inviability was one factor shaping the genetic composition of the F2 population in this inter-subgeneric hybrid. The conservation of microsatellite loci and synteny between Corymbia and Eucalyptus suggests there will be substantial value in exchanging information between the two groups.

How accurate are the marker orders in crop linkage maps generated from large marker datasets?

Marker ordering during linkage map construction is a critical component of QTL mapping research. In recent years, high-throughput genotyping methods have become widely used, and these methods may generate hundreds of markers for a single mapping population. This poses problems for linkage analysis software because the number of possible marker orders increases exponentially as the number of markers increases. In this paper, we tested the accuracy of linkage analyses on simulated recombinant inbred line data using the commonly used Map Manager QTX (Manly et al. 2001: Mammalian Genome 12, 930-932) software and RECORD (Van Os et al. 2005: Theoretical and Applied Genetics 112, 30-40). Accuracy was measured by calculating two scores: % correct marker positions, and a novel, weighted rank-based score derived from the sum of absolute values of true minus observed marker ranks divided by the total number of markers. The accuracy of maps generated using Map Manager QTX was considerably lower than those generated using RECORD. Differences in linkage maps were often observed when marker ordering was performed several times using the identical dataset. In order to test the effect of reducing marker numbers on the stability of marker order, we pruned marker datasets focusing on regions consisting of tightly linked clusters of markers, which included redundant markers. Marker pruning improved the accuracy and stability of linkage maps because a single unambiguous marker order was produced that was consistent across replications of analysis. Marker pruning was also applied to a real barley mapping population and QTL analysis was performed using different map versions produced by the different programs. While some QTLs were identified with both map versions, there were large differences in QTL mapping results. Differences included maximum LOD and R-2 values at QTL peaks and map positions, thus highlighting the importance of marker order for QTL mapping

Homologues of the maize rust resistance gene Rp1-D are genetically associated with a major rust resistance QTL in sorghum

As part of a comparative mapping study between sugarcane and sorghum, a sugarcane cDNA clone with homology to the maize Rp1-D rust resistance gene was mapped in sorghum. The cDNA probe hybridised to multiple loci, including one on sorghum linkage group (LG) E in a region where a major rust resistance QTL had been previously mapped. Partial sorghum Rp1-D homologues were isolated from genomic DNA of rust-resistant and -susceptible progeny selected from a sorghum mapping population. Sequencing of the Rp1-D homologues revealed five discrete sequence classes: three from resistant progeny and two from susceptible progeny. PCR primers specific to each sequence class were used to amplify products from the progeny and confirmed that the five sequence classes mapped to the same locus on LG E. Cluster analysis of these sorghum sequences and available sugarcane, maize and sorghum Rp1-D homologue sequences showed that the maize Rp1-D sequence and the partial sugarcane Rp1-D homologue were clustered with one of the sorghum resistant progeny sequence classes, while previously published sorghum Rp1-D homologue sequences clustered with the susceptible progeny sequence classes. Full-length sequence information was obtained for one member of a resistant progeny sequence class ( Rp1-SO) and compared with the maize Rp1-D sequence and a previously identified sorghum Rp1 homologue ( Rph1-2). There was considerable similarity between the two sorghum sequences and less similarity between the sorghum and maize sequences. These results suggest a conservation of function and gene sequence homology at the Rp1 loci of maize and sorghum and provide a basis for convenient PCR-based screening tools for putative rust resistance alleles in sorghum.

Functionally associated molecular genetic marker map construction in perennial ryegrass (Lolium perenne L.)

A molecular marker-based map of perennial ryegrass (Lolium perenne L.) has been constructed through the use of polymorphisms associated with expressed sequence tags (ESTs). A pair-cross between genotypes from a North African ecotype and the cultivar Aurora was used to generate a two-way pseudo-testcross population. A selection of 157 cDNAs assigned to eight different functional categories associated with agronomically important biological processes was used to detect polymorphic EST–RFLP loci in the F1(NA6 × AU6) population. A comprehensive set of EST–SSR markers was developed from the analysis of 14,767 unigenes, with 310 primer pairs showing efficient amplification and detecting 113 polymorphic loci. Two parental genetic maps were produced: the NA6 genetic map contains 88 EST–RFLP and 71 EST–SSR loci with a total map length of 963 cM, while the AU6 genetic map contains 67 EST–RFLP and 58 EST–SSR loci with a total map length of 757 cM. Bridging loci permitted the alignment of homologous chromosomes between the parental maps, and a sub-set of genomic DNA-derived SSRs was used to relate linkage groups to the perennial ryegrass reference map. Regions of segregation distortion were identified, in some instances in common with other perennial ryegrass maps. The EST-derived marker-based map provides the basis for in silico comparative genetic mapping, as well as the evaluation of co-location between QTLs and functionally associated genetic loci.

Isolation and characterization of novel microsatellite markers and their application for diversity assessment in cultivated groundnut (Arachis hypogaea)

Background: Cultivated peanut or groundnut (Arachis hypogaea L.) is the fourth most important oilseed crop in the world, grown mainly in tropical, subtropical and warm temperate climates. Due to its origin through a single and recent polyploidization event, followed by successive selection during breeding efforts, cultivated groundnut has a limited genetic background. In such species, microsatellite or simple sequence repeat (SSR) markers are very informative and useful for breeding applications. The low level of polymorphism in cultivated germplasm, however, warrants a need of larger number of polymorphic microsatellite markers for cultivated groundnut. Results: A microsatellite- enriched library was constructed from the genotype TMV2. Sequencing of 720 putative SSR-positive clones from a total of 3,072 provided 490 SSRs. 71.2% of these SSRs were perfect type, 13.1% were imperfect and 15.7% were compound. Among these SSRs, the GT/CA repeat motifs were the most common (37.6%) followed by GA/CT repeat motifs (25.9%). The primer pairs could be designed for a total of 170 SSRs and were optimized initially on two genotypes. 104 (61.2%) primer pairs yielded scorable amplicon and 46 (44.2%) primers showed polymorphism among 32 cultivated groundnut genotypes. The polymorphic SSR markers detected 2 to 5 alleles with an average of 2.44 per locus. The polymorphic information content (PIC) value for these markers varied from 0.12 to 0.75 with an average of 0.46. Based on 112 alleles obtained by 46 markers, a phenogram was constructed to understand the relationships among the 32 genotypes. Majority of the genotypes representing subspecies hypogaea were grouped together in one cluster, while the genotypes belonging to subspecies fastigiata were grouped mainly under two clusters. Conclusion. Newly developed set of 104 markers extends the repertoire of SSR markers for cultivated groundnut. These markers showed a good level of PIC value in cultivated germplasm and therefore would be very useful for germplasm analysis, linkage mapping, diversity studies and phylogenetic relationships in cultivated groundnut as well as related Arachis species.

Domestication to crop improvement: Genetic resources for Sorghum and Saccharum (Andropogoneae).

Background: Both sorghum (Sorghum bicolor) and sugarcane (Saccharum officinarum) are members of the Andropogoneae tribe in the Poaceae and are each other's closest relatives amongst cultivated plants. Both are relatively recent domesticates and comparatively little of the genetic potential of these taxa and their wild relatives has been captured by breeding programmes to date. This review assesses the genetic gains made by plant breeders since domestication and the progress in the characterization of genetic resources and their utilization in crop improvement for these two related species. Genetic Resources: The genome of sorghum has recently been sequenced providing a great boost to our knowledge of the evolution of grass genomes and the wealth of diversity within S. bicolor taxa. Molecular analysis of the Sorghum genus has identified close relatives of S. bicolor with novel traits, endosperm structure and composition that may be used to expand the cultivated gene pool. Mutant populations (including TILLING populations) provide a useful addition to genetic resources for this species. Sugarcane is a complex polyploid with a large and variable number of copies of each gene. The wild relatives of sugarcane represent a reservoir of genetic diversity for use in sugarcane improvement. Techniques for quantitative molecular analysis of gene or allele copy number in this genetically complex crop have been developed. SNP discovery and mapping in sugarcane has been advanced by the development of high-throughput techniques for ecoTILLING in sugarcane. Genetic linkage maps of the sugarcane genome are being improved for use in breeding selection. The improvement of both sorghum and sugarcane will be accelerated by the incorporation of more diverse germplasm into the domesticated gene pools using molecular tools and the improved knowledge of these genomes.

QTL analysis of ergot resistance in sorghum

Sorghum ergot, caused predominantly by Claviceps africana Frederickson, Mantle, de Milliano, is a significant threat to the sorghum industry worldwide. The objectives of this study were firstly, to identify molecular markers linked to ergot resistance and to two pollen traits, pollen quantity (PQ) and pollen viability (PV), and secondly, to assess the relationship between the two pollen traits and ergot resistance in sorghum. A genetic linkage map of sorghum RIL population R931945-2-2 x IS 8525 (resistance source) was constructed using 303 markers including 36 SSR, 117 AFLP™, 148 DArT™ and two morphological trait loci. Composite interval mapping identified nine, five, and four QTL linked to molecular markers for percentage ergot infection (PCERGOT), PQ and PV, respectively, at a LOD >2.0. Co-location/linkage of QTL were identified on four chromosomes while other QTL for the three traits mapped independently, indicating that both pollen and non pollen-based mechanisms of ergot resistance were operating in this sorghum population. Of the nine QTL identified for PCERGOT, five were identified using the overall data set while four were specific to the group data sets defined by temperature and humidity. QTL identified on SBI-02 and SBI-06 were further validated in additional populations. This is the first report of QTL associated with ergot resistance in sorghum. The markers reported herein could be used for marker-assisted selection for this important disease of sorghum.